On complex singularities of the 2D Euler equation at short times 



o 
o 

(N 



00 

Q 
U 



> 

cn 

cn 
o 

ON 

o 



W. Pauls' 

'Max Planck Institute for Dynamics and Self-Organization, Gottingen, Germany 
(Dated: March 18, 2009) 

We present a study of complex singularities of a two-parameter family of solutions for the two-dimensional 
Euler equation with periodic boundary conditions and initial conditions ?/;o( 21,22) = i^(p)cosp ■ z + 
i^(q) cos q ■ 2: in the short- time asymptotic regime. As has been shown numerically in W. Pauls et al, Physica 
D219, 40-59 (2006), the type of the singularities depends on the angle between the modes p and q. Here we 
show for the two particular cases of (j) going to zero and to n that the type of the singularities can be determined 
very accurately, being characterised by a = 5/2 and q = 3 respectively. In these two cases we are also able to 
determine the subdominant corrections. Furthermore, we find that the geometry of the singularities in these two 
cases is completely different, the singular manifold being located "over" different points in the real domain. 



I. INTRODUCTION 

One of the most important open problems in mathematical 
fluid dynamics is the question of well-posedness of the three- 
dimensional incompressible Navier-Stokes and Euler equa- 
tions (for the most recent reviews see |1, 2, 3]). Although 
the Navier-Stokes equation is more relevant from the physical 
point of view, certain phenomena related to the non-locality 
of the equations are partially obscured by the presence of the 
viscous term and are more conveniently studied in the invis- 
cid case. One of these phenomena is the depletion of non- 
linearity: generally flows behave much more tamely than one 
would expect based on simple dimensional arguments and 
tend to suppress singularities, approaching - at least locally - 
one of the numerous steady state solutions of the Euler equa- 
tion. 

Typically, inviscid incompressible flows observed in nu- 
merical simulations exhibit a much slower loss of regular- 
ity than that predicted by analytical estimates. However, it 
should be taken into account that numerical simulations aim- 
ing at studying the blow-up are mostly performed on periodic 
domains by employing the pseudo-spectral method with very 
smooth initial conditions, such as trigonometric polynomials. 
In contrast to numerical studies, analytical estimates aim at 
much larger classes of initial conditions for which the rate of 
loss of regularity may be much faster Indeed, in two dimen- 
sions Bahouri and Chemin [4] have constructed an example 
of a vortex patch initial condition with exponentially fast reg- 
ularity loss, thus proving that the corresponding estimate of 
the regularity loss is sharp. Therefore, it is important to clar- 
ify how the strength of depletion may depend on the character 
of the initial conditions. 

One possibility to assess the strength of depletion and its 
dependence on the initial conditions consists in analysing 
complex singularities' of solutions of the Euler equation: for 



real analytic initial conditions, solutions can be continued into 
the complex domain by extending the space variables to com- 
plex values. These solutions typicaly become singular at some 
complex locations constituting, according to the numerical ev- 
idence given in |6, 7], a smooth manifold. Depletion is espe- 
cially pronounced in the neighbourhood of the complex singu- 
larities: assuming that the singular manifold has co-dimension 
one, it follows from the incompressibility condition that non- 
linearity vanishes to the leading order in the direction normal 
to it. 

It has been found in [Isfl (in the following abbreviated as 
PMFB) by high-precision numerical simulations that this pro- 
liferation of depletion in the neighbourhood of singularities 
entails a significant dependence of the analytic properties of 
solutions on the initial conditions: it affects the geometry but 
also the type of the singularities. The singularities are non- 
universal: strong numerical evidence has been given in PMFB 
that in the neighbourhood of the singularities the vorticity be- 
haves like (here s is the distance to the singular manifold) 
with the exponent (3 depending on the choice of the initial con- 
dition. 

Let us now briefly describe some of the findings in PMFB. 
In that work, solutions coiTesponding to the initial conditions 
ipo{zi, Z2) — F{p) cos|) • z + F{q) cos q • z were considered 
in the short-time asymptotic regime. The exponent f3 was de- 
termined indirectly by analysing the asymptotic behaviour of 
the Fourier coefficients of the stream function along so-called 
diagonals (i.e. directions with rational slope) in the Fourier 
space at high wavenumbers. Asymptotically, the Fourier co- 
efficients were found to behave as product of an exponen- 
tial (which gives the leading order contribution) with an al- 
gebraic prefactor |k|~" (being the first subdominant contribu- 
tion). The exponent a is related to the rate of divergence of 
the vorticity (3 via a + f3 = 7/2 and can be determined nu- 
merically. The precision in determining the numerical value 



An additional reason for studying complex singularities is that possible 
emergence of real singularities in three dimensions can be detected by 
monitoring the dynamics of complex singularities 1 12], since, as has been 



shown by Bardos et al |^ and Benachour iTolfTTIl , a hypothetic real singu- 
larity can only arise when a singularity located in the complex domain hits 
the real domain. 
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of the exponent a was sufficient to insure that the type of the 
singularities changes depending on the initial conditions: two 
different initial conditions were found with a — 2.66 ± 0.01 
and a — 2.54 ± 0.01. It was also found that the numerical 
value of the algebraic prefactor exponent a depends on the 
angle between the vectors |) and q, but does not depend on 
their lengths. 

To achieve a better understanding of the non-universality 
of the type of complex singularities and of their dependence 
on the initial conditions a systematic study of the numerical 
values of a and their dependence on the angle between p and 
q is needed. However, such a study requires a more precise 
determination of the values of the exponent a as a function of 
the initial data, since for small values of the angle between |) 
and q the changes in the numerical value of a are of the same 
order as the errors in a. 

Unfortunately, for the cases which have been presented in 
PMFB a more precise determination of the exponent is not 
feasible with the resolution achieved there. The main rea- 
son for this is the presence of strong intermediate asymptotics 
and our lack of knowledge of the asymptotic expansion of the 
Fourier coefficients along diagonals in the Fourier space. 

Here we shall show that there exist initial conditions for 
which it is possible to achieve high precision for the numerical 
value of a and to determine the subdominant contributions: in 
the one case we find a = 5/2, in the other a = 3. In Sec- 
tion we introduced a rescaled form of the two-dimensional 
Euler equation, describe Taylor expansions of its solutions, 
their relation to complex singularities and to the short-time be- 
haviour of solutions. In Section|III]we give a detailed account 
of the short-time asymptotic regime, in particular with respect 
to the dependence on the angle between the vectors 43 and q. 
In Section IrV] we study the case of initial conditions with the 
angle between p and q going to zero. We show that the nu- 
merical value and the asymptotic expansion of the Fourier co- 
efficients can be determined very accurately. In Section [V] 
we study the limiting case of the vectors p and q becoming 
anti-parallel. We show that in contrary to the cases studied 
previously in PMFB the complex singularities are not located 
"above" a symmetry point. We also determine accurately the 
value of a and the asymptotic expansion of the Fourier coef- 
ficients. In Section|Vl]we make some remarks on the general 
case, when the angle between p and q is not restricted to any 
special value. Finally, we conclude with Section lVlIl 



11. RESCALED TWO-DIMENSIONAL EULER EQUATION 

A. Basic equation 

We begin by considering the two-dimensional Euler equa- 
tion 

atAV^ = J(^,AV'), (1) 

on the periodic domain — [0, 27r) x [0, 2-k). We chose the 
simplest periodic initial conditions with nontrivial dynamics. 



that is, trigonometric polynomials with only two harmonics 
^0(^1, ^2) = F{p)e-''^-^ + F(q)e-*''-^ + c.c. (2) 

Note that after a real translation we can assume that the initial 
amplitudes are real. Therefore the initial condition becomes 

V'o(^;i,z2) = 2/'(}3)cos|j-2 + 2F(q)cosq-2;. (3) 

Since we shall restrict ourselves to the initial conditions 
we introduce the auxiliary coordinates zi = p • 2;, Z2 = q • 2; 
and a rescaled time r = (|j A q)t. Then, writing the solution 
as 

^^(2:1,2:2;*) =V^(2i,z2;(pAq)t), (4) 
we obtain the Euler equation in the following form 

9rA(p,„)7/; = J(V', A(p_,)-0), (5) 
with initial condition 

^'0(^1,22) = 2/'(fi)coszi + 2F(q)cosz2, (6) 
where A(p q) denotes the Lapacian 

A(p,,) = |ppaJ + 2p.qaia2 + |qp92- (7) 

Note that in this formulation the vectors p and q do not nec- 
essarily have to be vectors in 1? , but can take any arbitrary 
values in M^. 

It is obvious that if ^"(^1 , ^2 ; i) is a solution of (HI with in- 
titial condition t/iq, then V'a('2i,'z2; t\) = A7/';^(zi, Z2; At) is 
the solution of the same equation with rescaled time t\ = Ar 
and the initial condition \tpQ. Thus, the only relevant parame- 
ters characterising the problem are the ratio A = F{q)/ F{p) 
of the two initial amplitudes F{p) and F{q) and the angle (p 
between the vectors p and q together with the ratio of lenghts 
77 — \q\/\p\ of these two vectors. The latter fact is immedi- 
ately seen by rewriting the Laplacian as 

A(p,q) = IpI A(0_^), (8) 

and obtaining 

5rA(^_^)l/' = J(V', A(^^^)l/'), (9) 

where 

1 2 2 

^ -di + '2 cos4>did2 + r]d2- (10) 



B. Perturbative approach 

From ( [Tol l it follows immediately that in the case — 1 
Equation (|9]) has a stationary solution. Since, as we have 
mentioned, the parameter 7/ can in principle take any values, 
for values of rj close to 1 we try a perturbative expansion in 
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£ = rj — 1. The leading order t/'q coincides with the initial 
condition 

i^o{zi,Z2) = 2 F{p) COS zi + 2F(q)cosz2, (11) 

whereas the next order V'l is the solution of the linearised Eu- 
ler equation with a source term 



2 F{p) F{q) sinzi sinz2- 



(12) 



The background flow ipQ is the well-known cat's eye flow. 
Note that neglecting the source term and the perturbation term 
J{ipi, ^ipipo) we obtain the shortened Euler equation which 
is just the passive scalar equation 



dtA4,ij^+v- V(A0Vi) ^ 0, 



(13) 



where the velocity field is determined by i) = 
{d2tpa, —di4>o)- Singularities of advected fields such as 
A^f/j]^ have been studied in [ 15, 16]. Since these singularities 
are determined by the singularities of the back-to-labels map 
and the flow is integrable, the singularities can be determined 
completely, see fis'l for the case A = 1 (F{p) ~ F{(\)) and 
lll6i] for the more general case < A < 1. 

Coming back to the linearised Euler equation we note that it 
has been extensively studied for various types of background 
flows. In particular, in |17] flows fo Kolmogorov's type have 
been studied (this type of flows, however, does not have any 
complex singularities). However, an analytic study of sin- 
gularities of solutions of the linearised time-dependent Euler 
equation has not been attempted here, but certainly merits a 
more detailed consideration. 



C. Taylor series expansion 

To calculate the solutions of Equation ^ we use the ap- 
proach based on Taylor series expansion [,21il : representing 
the solution as Taylor series in time 



(14) 



it is obvious that the functions -ip (z) satisfy the following 
recursion relation 

-r- -j(n+l) i Sr^-ffM"^) -r -r("-™)N .i^s 



m=0 



<1) 



For example, for ip we obtain 

2F{p)F{qy-^L_^cos{z, + Z2). 



(16) 



III general the functions 4' (zi , Z2) will have the form 

2m<n n— 2m 

^<"\zi,Z2)= E E 2{Fip)r^'~'--{Fiq)y-^^ 

m—0 cj— 1 

{Fj";,"^ cos((n + 1 - 2m - a) 21 +^22) + 
pin.jn) j,Qg^|-j^ + 1 — 2rn — a)zi — crz2) |-|- 



2m<n 



^ 2(F(p))"+'~'"(F(q))''"F,("^'") cos((n + 1 - 2m)z,] 



m—1 
2m<?i 



Y: 2(F(^.))-^'"(F(q))"+^-^'"Fi"-) cos((n + 1 - 2m)z2 



m—1 



(17) 



The coefficients fI"^''\ f1"'J^\ i^i^"'™' and F^''^"'^ depend 
only on the initial vectors |j and q, more precisely on the pa- 
rameters (f> and r/. For example, for t/'^^'' we have 



p(2,0) ^ 1 M 

+ 1 2 \p + q\^ 



\2p 



F 



(2.0) 



1 |q| 



F, 



2 \p- 

(2.0) _ _i |q| 



' \2p-q\ 
\p\^ \p + q\^- 



F_ 



:2,o) 



2 1^ 
1 |qP 



|P + 2q| 



\P\' \p-^\' 



F, 



(2.1) 



Fr 



(2.1) 



2 |^-q| 

l |qp-|yP 
2 IpP 

l |qp-|P 
2 |qP 



|p-2q|2 



(18) 



Inserting the representation (T% into Equation ( fTSl l we obtain 
a system of recursion relations for the coefficients. Since the 
resutling expressions turn out to be quite bulky we prefer not 
to list them here (see, however, [22] in a slightly different no- 
tation). 

Note that contrary to solutions of the original Euler equa- 

tion the higher-order terms 1/; (21,22) do not vanish trivially 
for (/) — > (or — > tt) if we take the limit keeping the ratio 
?7 of lengths of p and q fixed. However, it is not obvious that 
the limit (or </) — > tt) is well-defined. Indeed, if we 
try to evaluate expressions ( fTSl ) by putting directly, say = 
and 1] = 1 or 77 = 2, we obtain ill-defined expressions because 
of vanishing denominators. Thus, we have to be careful about 
specifying the limit 0-^0. Actually we are considering a 
double limit: the limit of the initial modes becoming parallel 
and the limit of long times, as can be seen from (|4|l. We shall 
come back to this question in Sectionlllll 

Let us finally make a small remark on the utility of the 
time Taylor series expansion: introducing the auxiliary co- 
ordinates ^1 = e~'^i and ^2 = e~*^^ we can write the so- 
lution in form of a triple Laurent series in r, ^1 and ^2- The 
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solution obtained in this form is only meaningful inside the 
domain of convergence of the corresponding series^. The do- 
main of convergence is determined by the complex singulari- 
ties of the function ^p{£,i,£,2', t): in particular, because of the 
regularity of two-dimensional flows we know that for real zi, 
02 there are no singularities for real values of r. However, 
there will typically be singularities for some complex-valued 
T*, which will prevent the series from convergence even for 
real-valued r's, when r > |r*|, see e.g. |23l1 . Thus, the power 
series representation is not well-suited for studying the tem- 
poral dynamics of the solution. Nevertheless, it lends itself 
very well for analysing space-time structure of the singulari- 
ties. Indeed, it is known that from the asymptotic behaviour of 
the coefficients of a series conclusions can be made concern- 
ing the geometry and the type of the complex singularities of 
V^(6,6;T),see|i|2i. 



D. Short-time asymptotics 

The calculation of the solution is greatly simphfied if we 
consider the asymptotic short-time behaviour as in IE S Htl • 
Indeed, at short times the Fourier coefficients corresponding 
to each cos(fciZi + k2Z2) can be represented as 

(19) 

for k2 > and 

(20) 

for k2 < 0. The wave vector ki can be assumed positive 
without any loss of generality. Then, in the asymptotic regime 
t ^ the solution will have the form 

?as(^i,^2) = f:E2(F(p))"+^-^(i^(q))Vx 

n=0 a=l 

cos((n + 1 - a)z, + az2)+ ^^^^ 
cos((n + 1 - a)zi — cfz2)^ 

Representing cos ( (n + 1 — a)zi + az2 and cos ( (n + 1 — a)zi — 
(JZ2 as a sum of two exponentials we can write the short-time 
asymptotic solution as 

V'as('2i7^2) = -^+(^1 + ilnr, Z2 + ilnr)+ 

1_ ^ (22) 

— ij}_{'zi + ilnr, Z2 — ilnr) + c.c, 



^ Unless we perform analytic continuation, which is the subject of an ongo- 
ing work. 



where 

^+{zi,Z2) = 
oo n 

n=0 (T=l 

(23) 

and, similarly 

i^-{zi,Z2) = 
OO n 

(24) 

As has been noted in the functions 

— ^p, (zi + i In T, '22 + i In r) , — {zi + ilnT,Z2 — i In t) 

T T 

are actually solutions of the Euler equation ^ corresponding 
to initial conditions 

F{p)e-''' + F{q)e-''^ (25) 

and F(|3)e"*^i + F(q)e'^^. Thus, at short times, the solution 
of the Euler equation (|4|i with initial condition ^ behaves as 
a linear superposition of four solutions with very simple two- 
mode initial conditions of the type ( |25] | 

Geometrically, the above observation can be explained as 
follows (see Fig. [T]): the singular manifold consists of four 
branches which at short-times are located very far away from 
each other Therefore, at short times every branch behaves as 
if it would not see the presence of other branches; the nonlin- 
ear interactions between separate branches are weak. 




FIG. 1: Schematic representation of the projection of the singular 
manifold onto the (j/i , j/2)-plane at short-times. Only the parts of the 
singular manifold closest to the real domain are sketched. The sin- 
gularities are represented in non-rescaled coordinates. Note that the 
vectors p and q determine the asymptotic directions of the various 
branches. 
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ni. SOLUTIONS OF THE EULER EQUATION IN THE 
SHORT-TIME ASYMPTOTIC REGIME 

A. Fourier space and pseudo-hydrodynamic formulations 

In this section we study short-time asymptotic solutions of 
the two-dimensional Euler equation (|4|. Since the coefficients 
i^"^"-* can be determined using the formula 



Fi";°)(p,q) = (-im'r(p,-q 



(26) 



we concentrate on the solution given by Here we assume 

that 7^ 0. In what follows we put for brevity f}^'^^ = 

The coefficients Gct""* can be calculated by means of the 
following recursion relation 



G(") = -i - V V 

n |(n+l-a)p + aq|2 /-'^^^ 
|((n-m) - {cr-T))p + ((T-T)qp 
[(m + l)(a - r) - (n - m)T]G(.™)G^'r™"'^ 

with initial values G^,"^ = 1 = g[°\ Itholds^ 



(27) 



(28) 



Similarly to what has been noted in PMFB06, the coefficients 
G^""* can be obtained explicitly 



G 



„1 '-\ap + q\^-\p\ 



(n) _ . -,x„ -L TT 
1 ^) 



(29) 



(7=0 

as well as the coefficients Gri"^ • 

Analogously to PMFB the short-time asymptotics can be 
formulated as a stationary hydrodynamical problem by intro- 
ducing the pseudo-hydrodynamic stream function 

oo n 

Gih^h) = e'-' + + ^ ^ Q{n)^{n+l^a)z^^az2^ (3Q) 
n—1 a—1 

In terms of the pseudo-hydrodynamic stream function 
G(zi , Z2) the solution of the Euler equation Q with the initial 
condition F{p)e~^'^^ + i^(q)e^'^2 can be written as 



- G[~izi + Ini + lnF(43), -iz2 + Int + InF(q)). 



(31) 



This pseudo-hydrodynamic stream function solves the time 
independent equation 



A(p,q)G = 

(9iG + 1) (92A(p.q)G) - {d^G - 1) (aiA(p,q)G). 



(32) 



' Easily seen by permuting k\ = n + 1 — a and k2 = u and using formula 



The initial condition is replaced by a boundary condition at 
infinity: for zi —00, Z2 — > —00 holds 

G{zi,h)^e~'' +e'\ (33) 

As has been noted in |8], Equation ( l32b can be written in a 
form which is closer to (01 by introducing an auxiliary func- 
tion 

iJ(zi, Z2) = G(zi, 22) + zi - Z2. (34) 
We then obtain the following pseudo-hydrodynamic equation 
A(p,,)i/- J(i7,A(p,,)i/). (35) 



B. Parameter-dependence of tlie solutions and the limiting 

cases (1)^0 and (ji ^ tt 

To get an idea about the properties of the solutions in 
dependence on the parameters we calculate the coefficients 
Ga^\4>, 77) symbolically as functions of 4> and up to ri = 8, 
the coefficient g'"^'' [cj), 77) being represented on Fig. |2] From 




FIG. 2: Coefficient G^P {<j>, rf) in dependence on and rj. 

the study of the low-order coefficients we can infer that: for (i) 
(f) close to the dependence of the coefficients on both (j) and 
77 is very smooth, whereas for (ii) </) close to tt the coefficients 
behave very irregularly. This suggests that the case (/> — > 
may have a reduced non-linearity in comparison to the case 
(j) TT. With this motivation in mind we are going to study 
the limits cj) ^ and (j) ^ tt. 

In the case </) ^ the short-time asymptotic solution is 
obtained in a straightforward fashion by putting (p ^ 0. The 
Fourier coefficients of the solution can be calculated directly 
from the recursion relation (fTSl l with = 0. Indeed, the re- 
cursion relation ( fTSl ) is well-defined because for (fi — the 
Laplacian {{n + 1 — cr) + y/V'^Y does not vanish. 

The case tt is much more delicate: the values of the 
coefficients depend on the procedure we choose to calculate 
them. Consider, for example, the behaviour of the coefficient 
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Cf^ ((?!), 7?) = Fl^f\p, q) given by As function of 

and 77 it is not continuous at 77 = 1, 77 = 2 and (j) — tt. Indeed, 
if we fix first the angle = tt we obtain the expression 



1 7] + 1 rj 



2 77- 1 ?7- 2' 



(36) 



which has poles^ at 77 = 1 and 77 = 2. On the other hand, fix- 
ing 77 = 1 or 77 = 2 we will obtain finite values by considering 
the limits 



lim G\ 

0— >7r 



(2)/ 



1) =0, 



limGf)(^,2) = -. (37) 



Thus, the numerical values of the coefficients for given {rj, (f>) 
depend on the limiting procedure ^ tt, 7/' — > r;. Since 
we are mostly interested in the dependence of the solutions 
on the angle cf), it seems natural to consider the limit of the 
coefficients obtained by fixing the ratio of lengths 77 first and 
then letting (p ^ tt. The actual implementation of this limit 
will be discussed in Section[V]and AppendixlB] 

Note that the ^ (or ^ tt) limit obtained by fixing 
77 can also be obtained for the two-dimensional Euler equa- 
tion ([5]) with real-valued initial conditions (|6]), see Appendix 
IB] The relation to the so-called hydrostatic limit [18^ ,19. , i20ll 
remains to be clarified. 



IV. CASE OF PARALLEL INITIAL MODES 

In this section we shall consider numerical solutions for the 
case (j) = 0. As has been found in PMFB, the value of the 
exponent of the algebraic prefactor does not depend on the ra- 
tio of the lengths of the initial vector rj. In the following for 
convenience we shall consider the case 77 = 2. We have cal- 
culated the numerical solution by using the recursion relation 
( fTSl l with a precision of about 70 significant digits and with the 
resolution iVmax = 2000. All numerical calculations in this 
and the two following sections have been done using multiple- 
precision calculation packages MPUN90 package ll24ll ) and 

MPFR m. 

We have studied asymptotic expansion of the coefficients 

G(")(cr) = G(fci, fca) (where /ci = 7^ + 1 - cr and fcs = cr) 
along half-lines 

\k\{cos6, sinO) — (fci, = 77i(p, q) 

in the Fourier space (later we shall refer to them as rational 
directions) for |fe — > cx). For this purpose we have used the 
asymptotic interpolation procedure of van der Hoeven ll26ll . 
see also |27]. In Appendix lAl we briefly remind the reader 
the main principles of numerical analysis of asymptotic ex- 
pansions. 

We have chosen rational directions such that p,q < 10, 
where p and q are mutually prime numbers. This gives 55 



' Note that if we fix a rational rj, then there exist n and cr such that the Fourier 
coefficient F^"^ (tt, ri) has a pole at rj. 



different rational directions. First, after six stages of interpo- 
lation we have identified the leading order term and the two 
successive sub-leading order terms 



G{k,9) - G(6')fc-"('^)e-'^('')^ 



(38) 



Fig. [3] shows the discrepancy of the scaling exponent a, de- 
termined from the interpolation value at the sixth stage, from 
the value 5/2 as a function of the angle 9. It is seen that the 
numerically determined Ofnum differs from 5/2 by less than 
10^^. We remind the reader that a theoretical argument has 
been given in PMFB suggesting that the true value of a does 
not depend on 9. 













' 6/5 - a^'"" + 
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FIG. 3: Difference between the conjectured extrapolated value 
3/a — 6/5 for a — 5/2 and the numerical value 3/anum obtained 
at the sixth stage of the asymptotic interpolation procedure. 

To confirm that the algebraic prefactor exponent a = 5/2 at 
the next step we have tried to identify higher-order corrections 
to ( [38l l. Indeed, if the higher-order corrections are known, 
we will be able to estimate the parameters S{9), a and C{9) 
more precisely. Furthermore, to obtain some information on 
the analytical structure of the solutions it is advantageous to 
know their asymptotic expansion at high orders. 

We remark that at the sixth stage of the interpolation proce- 
dure the numerical data are of the form 



g^^-* (rn) = h uj{m), 

a 



(39) 



where 777, is the number of the point m{p, q) on the line with 
rational direction (p, q). Continuing the procedure beyond the 
sixth stage we have first identified the nature of the remainder 
ijj{m) by the procedure described in the AppendixlAl 

The identification has been performed for the coefficients 
corresponding to the most populated direction {p,q) = (1,1), 
which is the direction with the most pronounced asymptotic 
behaviour We have found that the correction is of the form 
^(777) ^ \/ni?. This correction corresponds to the following 
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form of the asymptotic expansion of the Fourier coefficients 

G{k,e) = 
C(e')fc-^e-'^W'- 



ai(g) a2(g)lnfc fl 

k P +^U2 



(40) 



We note that the same asymptotic expansion is found for solu- 
tions of the Hnearised Euler equation considered in the short- 
time regime, see AppendixIC] 

It is very unUkely that the functional form of the asymptotic 
expansion depends on the direction (p, q) in the Fourier space. 
However, for large p and q the number of points on the half- 
Unes with the rational direction (p, q) is not very high, and 
therefore the asymptoic behaviour may not be attained. There- 
fore the values of the coefficient 02 {9} that we have measured 
using the asymptotic interpolation at thirteenth stage exhibit a 
quite large scatter Due to this, for determining S{9) and C{6) 
(represented on Figs.|4]and|5]) we have taken into account only 
the correction term ai {6) / k. 



FIG. 4: Leading-order asymptotics: 5(9) in lin-log coordinates. 

As we have stated before, averaged over all direction the 
absolute error in determining a is of the order 10~^. However, 
restricting ourselves to the most populated rational directions, 
we can vastly increase the precision in determination of a. 
One possibility is for example to use one of the convergence 
acceleration algorithms such as the p-algorithm for the data at 
sixth stage ( [39] ). The choice of the algorithm is fixed through 
the form of the correction term being proportional to l/m?, 
see ll28ll . For example, for the most populated direction (1,1) 
the acceleration algorithm stabilises at 5/2 with an error of the 
order 10^^^ after 50 iterations. 



V. CASE OF ANTI-PARALLEL INITL\L MODES 

In this Section we study numerically the case of antiparallel 
modes ^ tt: as we have seen before, this limit cannot be 
analysed directly by taking (f) — tt. Instead, we fix 77, taking 
the limit cf) ^ tt of the recursion relation (fTSl l. 



C(e) 



"+++ + 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



FIG. 5: Second subleading term C(S). 



Here we shall only consider the case^ rj = 2: the recursion 
relation is undefined for n + 1 — a = 2a which we have to 
analyse separately. It turns out that both the sum on the right- 
hand side of ( fTSl l and the Laplacian ((n+ 1 — cr) — 2a)^ vanish 
in such a way that their ratio yields a finite limit. The detailed 
calculation is not difficult but somewhat cumbersome; its de- 
tails are given in Appendix IB] 

Numerical calculation of the coefficients reveals that the 
properties of the (f) 7r-limiting solution are very different 
from the properties of the solutions which have been observed 
before, e.g. from those of </> = or </> = 7r/2. Remember that 
in these cases the Fourier coefficients had alternating sings 
which means that the part of the singular manifold being clos- 
est to the real domain is located "above" the point (tt, 0) in the 
real space. 

In the present case the situation is more complicated. 
Firstly, it turns out that the Fourier coefficients G^"-' (cr) van- 
ish for (7 < n/3 and thus the support of the solution in 
the Fourier space is bounded by the half-lines fci = 1 and 
= (fci — l)/2. Secondly, as can be seen on Fig. |6l in 
the angular sector approximately between k2 — (fci — l)/2 
and /c2 ~ 3fci the coefficients can have both positive and neg- 
ative signs which are distributed in an interference pattern. 
And thirdly, we find that in the sector roughly between the 
half-lines fc2 = 3 fci and fci 1 the coefficients are strictly 
positive. This means that one part of the singular manifold is 
located over the point (0, 0), whereas the other part is found 
above some curve in the real space. 

The geometry of the singular manifold can be studied by 
a two-dimensional generalisation of the BPH-method, dis- 
cussed in 1 22], which is, however, out of the scope of the 
present article. Note that to study the type of the singularities 
it is sufficient to consider the sector with the positive coeffi- 
cients, applying to them the asymptotic interpolation proce- 



^ It is not difficult to extend the calculation to other values of »?, for example 
T] = 3 01 r] = 4. 
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FIG. 6: Distribution of the signs of the Fourier coefficients in the 
in the limiting case <j) n. All Fourier coefficients vanish for a < 
n/3. 



dure used previously. 

Here again for rational directions {p,q), q/p > 3 the lead- 
ing order term and the two successive sub-leading order terms 
are 



(41) 



and the data at the sixth stage of the asymptotic interpolation 
procedure have the form 



(to) = h uj^m). 



(42) 



Continuing the interpolation procedure up to the thirteenth 
stage we found that for slopes q/p bigger than 4 the correction 
can be determined very accurately and it is of the form 1 /m*. 
The asymptotic expansion is just an expansion in l/|fe| 



G{k,9) = 



-s{e)k 



k k"^ k^ 



(43) 



To increase the precision in the determination of a we use the 
p-algorithm to accelerate the convergence. We find that the 
constant 3/a at sixth stage stabilises at the value 1 with more 
than 30 digits. 

The value a = 3 means that the vorticity diverges as s~^/^ 
in the neighbourhood of the singular manifold. Note that for 
the two-dimensional Burgers equation with initial conditions 
that are given by the velocity field determined from the stream 
function (O, the derivative of velocity field diverges at the 
same rate. Also, the form of the asymptotic expansion ( |43] | 
(expansion in inverse powers of k) is identical to the asymp- 
totic expansion we would obtain for the Burgers equation. 



Thus, we conjecture that the complex singularities of the so- 
lutions of the Euler equation in the case <j) ^ tt are of the 
same type as the corresponding singularities of the Burgers 
equation. 



VI. SOME REMARKS ON THE GENERAL CASE OF 

< cjXn 

We have seen that the geometry as well as the nature of sin- 
gularities are very different in the two limiting cases tt — > 
and tt: in the former case the singularities closest to 
the real domain are located over the symmetry point (vr, 0) 
whereas in the latter case they lie over the point (0, 0) and a 
certain curve in the real domain. Furthermore, in these two 
cases the values of the scaling exponents (a = 5/2 in the for- 
mer and a = 3 in the latter case) and the functional forms of 
the asymptotic expansions are different. Therefore, it is not 
surprising to find that the general case < (p < ir exhibits 
strong intermediate asymptotics, making it difficult to deter- 
mine the actual high-fc behaviour of the Fourier coefficients. 
In particular this impedes a precise determination of the scal- 
ing exponent a and of the functional form of the asymptotic 
expansion. 

The easiest way to see the interplay between the 
and (f) ^ TT limiting behaviour is to look at the case when 
is very close to tt. First of all, in this case the Fourier co- 
efficients have non-vanishing values almost everywhere. An 
inspection of the sign distribution of the Fourier coefficients 
(Fig-13 shows the existence of three different regions: the first 




90 100 



FIG. 7: Sign distribution for an angle close to tt, here — 24/257r. 

region is roughly the one between the ki axis and the half-line 
with the slope 1/2. The Fourier coefficients in this sector have 
alternating signs. The corresponding part of the singular man- 
ifold is therefore located over the point (tt, 0). Note that in the 
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case = TT the Fourier coefficients in this sector vanish. The 
second region corresponds to the sector with the interference 
pattern in the case = tt. In the third region between the half- 
Hne with the slope 3 and the k2-axis the Fourier coefficients 
are strictly positive. The corresponding part of the singular 
manifold is therefore located over the point (0,0). The sin- 
gular manifold has a quite interesting geometry, however, its 
detailed study is out of the scope of the present article and will 
be presented elsewhere. 

When the angle (f> is decreased from tt to the region with 
alternating sings proliferates until all the Fourier coefficients 
become alternating in their signs. The threshold value of the 
angle is conjectured to be^ (j) — arccos(— 1/4). 

In view of the complicated geometry of the singular man- 
ifold it is not surprising that a precise determination of the 
scaling exponent is very difficult and needs very high res- 
olutions. One of the main factors reducing the precision in 
the determination of a{(f)) are the subdominant contributions 
which strongly depend on the directions in the Fourier space 
in which the Fourier coefficients are evaluated^. To increase 
the precision one would have to take into account these sub- 
dominant corrections. 

However, with resolutions of order A^max — 2000 we 
have not succeeded in determining the functional form of the 
asymptotic expansion of the Fourier coefficients and there 
could not obtain the numerical values of the scaling expo- 
nents with a precision comparable to the one found in the 
cases 0^0 and (f> ^ tt. Therefore, a systematic study of 
the dependence of the scaling exponents on the angle is not 
feasible at the moment. 



VII. CONCLUSIONS 

In the present paper we have shown that, even for the sim- 
plest initial conditions with non-trivial dynamics, complex 
singularities of the solutions of the two-dimensional Euler 
equation can be of very different nature, corresponding to the 
values of the scaling exponent a = 5/2 and a = 3. 

This has been achieved by identifying two particularly sim- 
ple cases for which the exponents and therefore the type of 
the singularities can be determined very accurately. Further- 
more, we were able to determine the subdominant terms. This 
increases the certainty in the numerical value of the scaling 
exponents far beyond the cases which have been analysed in 
PMFB. 

In the limiting case — > it has been found that the value 
of the scaling exponent a coincides with the value of the scal- 



* This threshold value applies only to the case rj = 2. In the more general 
case its value is arccos(~l/2»)). The thi'eshold values of have been 
estimated using symbolically calculated expressions for the Fourier coeffi- 
cients. 

^ This introduces an artificial dependence of the scaling coefficients obtained 
numerically on the direction 8 in the Fourier space. As we have mentioned 
before, it can be shown analytically that the value of the scaling exponent 
does not depend on 8. 



ing exponent for the solutions of the linearised Euler equa- 
tion. Furthermore, we have numerical evidence that the struc- 
ture of the asymptotic expansion in the Fourier space is the 
same in both cases. Thus, non-linearity is suppressed in this 
case, which constitutes an example of a dynamically nontriv- 
ial completely depleted flow. 

Actually, the case 0^0 can be treated by means of asymp- 
totic expansions near the edges, i.e. the regions for which 9 is 
either close to or close to tt. In particular, the value a ^ 5/2 
can be found both by Maple assisted calculations and asymp- 
totic arguments^. 

In the case (/> — > vr the scaling exponent a = 3 is the same as 
for solutions of the two-dimensional Burgers equation. Also 
in this case the structure of the asymptotic expansion coin- 
cides in both cases. This case is an example of a flow with 
a very strong nonlinearity, probably the strongest possible for 
the initial conditions of trigonometric polynomial type. 

The difference between the two cases cj) ^ and tt 
is reflected not only in the type of the singularities: we have 
found that in these two cases the geometry of complex singu- 
larities is completely different. In the case cf) all (except 
one) Fourier coefficients of the solution have regularly alter- 
nating signs which in the physical space is expressed by the 
fact that the part of the singular manifold closest to the real 
domain is located "over" the symmetry point (7r,0). In the 
case </) ^ TT only a part of the Fourier coefficients has regular 
signs being positive. The remaining coefficients do not have 
such regularity in the signs. The physical space counterpart 
means that the part of the singular manifold closest to the real 
domain is located partially "over" the point (0, 0) and partially 
over some curve in the real domain. 

For intermediate angles cf) lying between and vr, for exam- 
ple in the SOC case = 7r/2 we have not been able to de- 
termine the subdominant terms in the asymptotic expansion. 
Therefore, the uncertainty in the exact value of the scaling ex- 
ponent, in particular in dependence on the angle </> is quite 
high. 

What do our findings mean in consideration of the question 
of well-posedeness of the three-dimensional Euler equation? 
Firstly we note that the short-time asymptotic analysis of the 
kind presented here is easily extendable to three dimensions. 
Some preliminary results in this direction have been obtain 
in 1 22] confirming non-universality of the nature of complex 
singularities in three dimensions as well. Secondly one very 
important unresolved issue remains: does the nature of com- 
plex singularities change when we go beyond the short-time 
asymptotics? Although some indications have been given in 
PMFB that the nature of singularities does not change we have 
no definitive answer yet. 

We expect that the extent to which the non-linearity is de- 
pleted is strongly dependent on the initial conditions which 
we can also view as a dependence on the large-scale struc- 
tures in the flow. We cannot exclude that even if solutions of 
the three-dimensional Euler equation generically stay regular. 



A. Gilbert, private communication. 
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there exist initial conditions - such as the Kida-Pelz - with 
very strong non-linearity which could produce a finite-time 
blow up. 

Finally we note that the algebraic approach of the type 
used here can be extended to the viscous case. The short- 
time asymptotic regime of the inviscid case corresponds in 
the Navier-Stokes case to the asymptotic regime of small 
Reynolds numbers. 
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g{m)g{m — 2)/{g{m — 1))^ and D taking the difference 
g{m) g{iTi) — g{m — 1). After applying these sequence 
transformations the leading order of the obtained sequence to- 
gether with corrections will be 3/a + cj (m) . 

The main problem in studying the asymptotics of the 
Fourier coeffcients is the correct identification of the correc- 
tion term LLj{m). Indeed, once its form is known we can draw 
conclusions about the further sub-leading order terms in the 
asymptotic expansion iA.2\ of g{m). If Lo(m) is to the lead- 
ing order of the type m^^ , as in the case of the limits 0-^0 
and TT, we can identify the constant 7 by applying to the 
data after (IA.3I I the following transformations 

D, I, R, D, I, D, D. (A.4) 

The leading order term after these transformations will have 
the form 2/(7 + 1). In the cases and — > vr the 

constant 7 can be identified as 7 = 2 and 7 = 4 respectively. 



APPENDIX A: ASYMPTOTIC EXPANSIONS OF 
FOURIER/TAYLOR COEFFICIENTS 

Let us consider a two-dimensional array G{ki,k2) with 
ki,k2 > 0. One possibility to study the asymptotic behaviour 
of G{ki,k2) at infinity is to proceed by simple reduction, 
analysing the behaviour of the one-dimensional arrays^ 



5(p,g)("^) = G{mp,mq), 



(A.l) 



for relatively prime positive numbers p, q. An asymptotic 
expansion of gi^p^q){m) for m 00 yields then an asymp- 
totic expansion of G{ki, k2) in polar coordinates (/ci, /c2) — 
I k I (cos 0, sin 6) for | fc | . Note that other approaches have been 
suggested in |26]. They are not discussed here. 

The asymptotic expansion of g(p,q) {m) is found by apply- 
ing the asymptotic interpolation procedure 1^2^. Due to the 
presence of complex singularities the leading order and the 
two first sub-leading order contributions are 



(A.2) 



see also l8l|22t]. 

The main idea of the asymptotic interpolation procedure 
is to find for .9(p,g) such sequence transformations - stages 
of the interpolation procedure - which would "strip off" the 
sequence g(p.q) the leading order terms e^'^tpi)™ and m^" 
without knowing the precise value of the constants 5(^p.q) and 
a. One possible choice of transformations is the one which 
has been given in iIztIi : 



SR, -D, I, D, D, D, 



(A.3) 



where the sequence transformations are: I taking the inverse 
g{ra) l/g{m), SR taking the second ratio g{m) 



' The array g{mp,mq) is closely connected to the so-called diagonal of a 
double Taylor series, see (5^. 



APPENDIX B: PARALLEL AND ANTI-PARALLEL MODES 
LIMITS OF THE TWO-DIMENSIONAL EULER EQUATION 



1. Two-dimensional Euler equation in the < 



limit 



We begin by considering the two-dimensional Euler equa- 
tion (|9]l with initial conditions (|6]l. Putting formally = (or 
(j) = TT)we obtain 



1 — — 2— 

dt{—di±y^d2) V'form 



J 



(V'form, di ± 07 92) Vform) > 



(B.l) 



where the initial condition remain unchanged and the signs 
+ and — corresponds to the cases cj) — Q and = tt. Un- 
fortunately this equation is ill-posed, because, contrarily to 
the Laplacian A(0 ^) in the case < < tt, the operator 
((1/y^) di ± y/rjd2)'^ has a non-vanishing kernel. 

Here we show that if we first fix r/ and then let (/> — > 
we obtain a well-defined equation for the 0-^0 limit of (|9]l. 
Consider the two-dimensional Euler equation in the Fourier 
space 



dt'ip[ki,k2) 
1 



\pki + q/c2| 



l+l'=k 



(B.2) 



Clearly, for p = (1,0) and q = (ry,0), which corresponds 
to the limit 0^0 the denominator on the right-hand side 
vanishes. However, as we shall see now, the numerator also 
vanishes in such a way that the ratio of the both gives a finite 
number in the limit 0^0. 

We rearrange the sum in such a way that the term 
'0(/i, l2)i'{l'i, I2) appears only once for every I. Thus, we ob- 
tain terms of the form 

I A i {\pl[ + q/^|2 - \ph + qhl^) i^ihMm'iA)^ (B-3) 
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where I + I' = k. Noting that 
\pl[ + ql'2\^~\ph + ql2f = 

1 \2 /I \ 

-^ki + JTik2 - 2fcifc2(l - COS0) - 2li[-ki + fcs + 
^ / / 

2ZiA;2(l - cos(/)) - 2l2{ki + 7]k2) + 2l2ki{l - cos(p), 



and 



\pki + qk2\ 



(B.4) 



-^ki + y/nk2Y -2kik2{l-cos(l)), (B.5) 



we see that when ki + rjk2 = 0, it follows 

\pl[ + q^JiP - \ph + q^2p A:iA:2 - ^ife - ^2^1 



l^j/ci + qk2 



kika 



^0, 
(B.6) 

for fci, A;2 7^ 0. The only case where the above remark does 
not apply is when there exists an I such that 21 — k, in which 
case we obtain 



Ipi'i + 

\pki + qfcal 



'2? 



(B.7) 



Thus, the expression on the right-hand side of jB.lj is well- 
defined in the limit <j> Q. The analytic properties of this 
limiting solution and relation to solutions of the so-called hy- 
drostatic equation certainly merit a further study. We also note 
that a similar limiting solution can be constructed for the two- 
dimensional Navier-Stokes equation. 



2. Short-time recursion relation in the < 



TV limit 



and has thus a finite value. 

Indeed, assume that rii, n[, n2 and n^, are integers such that 
ni + n[ = 2a and 7^2 -|- rij = o. Then it holds 

|niJJ + n2qp - Irt'ip + nJjiP = 4(nin2 - n'^n2)(l + cos (/)), 

(B.ll) 

where we have used the expressions p = (1,0) and q — 
2(cos (p, sin (f). Similarly we obtain 



(n + 1 - cr)p + crq|2 = 8a^{l + coscj)), (B.12) 



and therefore the term 1 + cos 6 is cancelled. 



APPENDIX C: PERTURBATIVE STUDY OF THE 
SHORT-TIME ASYMPTOTICS 



Here we present a simple linear model which turns out to 
have the same asymptotic structure as the solution of the short- 
time Euler equation in the limit = 0. 

As we have noted in Section HTB] for small values of e = 
77 — 1 we can apply a perturbative ansatz to Equation ( l32b . 
which will give us to the first sub-leading order 

(1 - e^^)aiA(^,i)Gi + (1 + e'')d2\^.i)Gi ~ A(^,i)Gi + 
e^^diGi)^c^'id2Gi) = -2e''e'\ 

(C.l) 



Let us consider the limiting form of the recursion relation 
( fTSl l in the case 77 = 2 and (p ^ tt. Firstly, we remark that 
in the sum on the right-hand-side of this recursion relation the 
terms containing the product Gr^'G^'^L,™ appear twice: 
firstly for {m, r) and secondly for {n — m — I, a — t). Thus 
we can regroup the recursion relation in such a way as that the 
product Gr'"''G^"Lj.™^'' appears in only one summand, which 
has the form 

[(m + l){a -T)-{n- m)r\ G^r^G^:ir'\ ^^'^^ 
where 

K(^^^){n,a\m,T) ^ \ {(n ~ m) ~ (a ~ t))p + (cr -r)qp- 
| (m + 1 - r)p + rqp 

(B.9) 

Secondly, since for -q — 2 the recursion relation ( fT5] ) be- 
comes undefined only when cr = 2cr, we concentrate on 
the (f) — > TT-limit of the recursion relation only for n — ia —1. 
As we shall see now, when (j) ^ tt 

1 , , , 1 

■7-, ^ 7TT K(o a\in,a\m,T) = — - 

|(n + l-a)p + aq|2 (P^")^ ' I ' ^ 2a2 ^^^^^ 

I ((n — 77^) — (cr — t)) (ct — t) — (m + 1 — t)t^, 



Neglecting the source term and the terms e^^(9iGi) — 
e^^ (02 Gi) we obtain the passive scalar model which has been 
introduced in f^. For this model the singular manifold can 
be described analytically, and furthermore, we can determine 
analytically the type of singularities of the field AiGi, which 
turn out to be poles. So far we have not been able to handle 
Equation ( IC.ll ) satisfactorily by analytic tools. Therefore we 
have studied solutions of this equation by using high-precision 
numerics (for these calculations we have used a precision of 
about 70 digits and a resolution iVmax ~ 3000.). 

We have found that the presence (or absence) of the source 
term and of the term e^^(9iGi) — e^'^{d2Gi) changes nei- 
ther the location of the singularities nor the value of algebraic 
prefactor exponent a. Using the asymptotic interpolation pro- 
cedure of we were able to determine the asymptotic expansion 
of the solution in the Fourier space, which has the same form 
as the asymptotic expansion ( l40l i of the short-time Euler equa- 
tion in the limit cj) ^ 0. 

Note that the form of singular manifold and the scaling 
exponent 5/2 (or 1 in the physical space) do not depend on 
the parameter (j>. Thus, the perturbative expansion predicts a 
universal exponent and thus, gives wrong predictions when 
<f> ^ 0- Still, in the case = it gives the correct form of the 
asymptotic expansion. 
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